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Abstract 

This paper considers the use of singular perturbation approximations for a class of linear quan- 



[■"-«. ' turn systems arising in the area of linear quantum optics. The paper presents results on the physical 

(N '. 

reaUzability properties of the approximate system arising from singular perturbation model reduction. 



I. Introduction 

The modelling and control of quantum linear systems is an important emerging application 
area which is motivated by the fact that quantum mechanical features emerge as the systems 
being controlled approach sub-nanometer scales and as the required levels of accuracy in control 
and estimation approach quantum noise limits. In recent years, there has been considerable 
C^^ I interest in the feedback control and modeling of linear quantum systems; e.g., see ifTI- lfTTl . 
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Such linear quantum systems commonly arise in the area of quantum optics; e.g., see [USII - lEOll . 
The feedback control of quantum optical systems has applications in areas such as quantum 
/\f • communications, quantum teleportation, and gravity wave detection. In particular, the papers 

H ; 

— - ■ [HI, [[T5l - [[T7l have been concerned with a class of linear quantum systems in which the system 

can be defined in terms of a set of linear quantum stochastic differential equations (QSDEs) 
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expressed purely in terms of annihilation operators. Such linear quantum systems correspond to 
optical systems made up of passive optical components such as optical cavities, beam-splitters, 
and phase shifters. The main results of this paper apply to this class of linear quantum systems 
for the square case in which the number of outputs is equal to the number of inputs. 

This paper is concerned with the use of singular perturbation approximations in order to obtain 
reduced dimension models for the class of linear quantum systems under consideration. Singular 
perturbation approximations are is widely used for obtaining reduced dimension models for 
classical systems; e.g., see fTH . In the case of quantum systems, a reduced dimension model may 
be desired for a quantum plant to be controlled in order to simplify the controller design process 
which can be very complicated using existing quantum controller design methods such as the 
quantum LQG method of ^. Another application of model reduction for linear quantum systems 
arises in the case of controller reduction where a reduced dimension controller is obtained from 
a high order synthesized controller. In the case of coherent quantum control such as considered 
in jSl, (SI, [fTTl . the controller is required to be a quantum system itself and thus the reduced 
dimension system must be physically realizable. 

In the physics literature, a commonly used technique in the modeling of quantum systems 
is the method of adiabatic elimination, which is closely connected to the singular perturbation 
method in linear systems theory; e.g, see [|22ll - [|25l . The papers ll22ll - [|25l also consider the issue 
of convergence of these singular perturbation approximations. In this paper, we consider the 
properties of the singular perturbation approximation to a linear quantum system from a linear 
systems point of view; e.g., see ll26l for a detailed description of singular perturbation methods in 
linear systems theory including error characterization in both the time and frequency domains. In 
particular, we are concerned with the physical realizability properties of the singular perturbation 
approximation to a linear quantum system. The issue of physical realizability for linear quantum 
systems was considered in the papers [|5]|, JU, lfT6l . ifTTl . This notion relates to whether a given 
QSDE model represents a physical quantum system which obeys the laws of quantum mechanics. 
In particular, the results of the papers (51-113, lfT6l show that the notion of physical realizability 
enables a direct connection between results in quantum linear systems theory and linear systems 
theory. In applying singular perturbation methods to obtain approximate models of quantum 
systems, it is important that model obtained is a physically realizable quantum system so that it 
retains the essential features of a quantum system. Also, if the approximate model of a quantum 



plant is to be used for controller synthesis, the controller synthesis procedure may need to exploit 
the physical realizability of the plant model. In addition, if model order reduction is applied to a 
coherent feedback controller which is to be implemented as a quantum system, then this reduced 
order controller model must be physically realizable. 

In the paper [|T6ll . the notion of physical realizability is shown to be equivalent to the lossless 
bounded real property for the class of square linear quantum systems under consideration. This 
property requires that the system matrix is Hurwitz and that the system transfer function is unitary 
for all frequencies. The main result of this paper shows that if a singularly perturbed linear 
quantum system is physically realizable for all values of the singular perturbation parameter, 
then the corresponding reduced dimension approximate system has the property that all of its 
poles are in the closed left half of the complex plane and its transfer function is unitary for all 
frequencies. These properties indicate that in all but pathological cases, the singular perturbation 
approximation method will yield a physically realizable reduced dimension system. In addition, 
an example is given showing one such pathological system in which the singular perturbation 
approximation is not strictly Hurwitz. 

The paper also presents a result for a special case of the singularly perturbed linear quantum 
systems considered in this paper. This special case corresponds to singular perturbations which 
arise physically from a perturbation in the system Hamiltonian. In this case, the result shows 
that the corresponding reduced dimension approximate system is always physically realizable. 
This result can in fact be derived from the nonlinear quantum system results presented in the 
papers [|23l . [|25l . However, we have included this result, along with a straightforward proof, for 
the sake of completeness. We have also included an example of a singularly perturbed linear 
quantum optical system which fits into the subclass of singularly perturbed quantum systems 
for which this result applies. This example illustrates how such singularly perturbed quantum 
systems can arise naturally in physical quantum optical systems. 

The remainder of this paper proceeds as follows. In Section |Ill we define the class of 
linear quantum systems under consideration and recall some preliminary results on the physical 
realizability of such systems. In Section Hill we consider the singular perturbation approximation 
to a linear quantum system. We first present a result for the class of singularly perturbed linear 
quantum systems under consideration which relates to the lossless bounded real property. We then 
consider a special class of singular perturbations which is related to corresponding perturbations 



of the quantum system coupling operator and Hamiltonian operator. We present a result which 
relates to this class of singular perturbations and shows that the corresponding approximate 
reduced dimension system is guaranteed to be physically realizable. In Section |IVl we present 
a simple example from the field of quantum optics to illustrate the proposed theory. In Section 
rvl we present some conclusions. 

II. A Class of Linear Quantum Systems 

We consider a class of linear quantum systems described in terms of the annihilation operator 
by the following quantum stochastic differential equations (QSDEs): 

da{t) = Fa{t)dt + Gdu{t); 

dy(t) = Ha{t)dt + Kdu{t) (1) 

where F E C"^", G E C"^™, H E C'"^" and K E C™^™; e.g., see 0, Ml, El, lEl, [l20ll. 
Here a(t) = [ai(t) ■ ■ ■ a„(t)] is a vector of (linear combinations of) annihilation operators. The 
vector u{t) represents the input signals and is assumed to admit the decomposition: 

du{t) = (3u{t)dt + du{t) 

where u{t) is the noise part of u(t) and /^^(t) is an adapted process (see [27] and [jlSlD . The 
noise u{t) is a vector of quantum noises. The noise processes can be represented as operators on 
an appropriate Fock space (for more details see |l27|). The process /3„(t) represents variables of 
other systems which may be passed to the system © via an interaction. More details concerning 
this class of quantum systems can be found in the references lfT6ll . ^. 

Definition 1: (See [|7]|, lfT6l . ifTTl .) A linear quantum system of the form ([T]) is said to be 
physically realizable if there exists a commutation matrix = 0^^ > 0, a coupling matrix A, a 
Hamiltonian matrix M = M\ and a scattering matrix S such that 

F = -e(iM + -A^A\; 

a = -QA^S; 

H = A; 

K = S (2) 

and S^S = L 



Here, the notation "^ represents conjugate transpose. In this definition, if the system ([T]) is 
physically realizable, then the matrices S, M and A define a open harmonic oscillator with 
scattering matrix S, coupling operator L = Aa and a Hamiltonian operator Ti = a) Ma\ e.g., see 
[Ell, [|27l, [l29l and [El. This definition is an extension of the definition given in [[5l. [[T6l. l[T7l 
to allow for a general scattering matrix S\ e.g., see [Q. 

The following theorem is a straightforward extension of Theorem 5.1 of [[T6l to allow for a 
general scattering matrix S. 

Theorem 1: (See [16].) A linear quantum system of the form Q) is physically realizable if 
and only if there exists a matrix 6 = 6^^ > such that 

FQ + QF^ + GG^ = 0; 

G = -QH^K; 
K^K = L (3) 

In this case, the corresponding Hamiltonian matrix M is given by 

M = ^ (e~^F - F^e^i) , (4) 

the corresponding coupling matrix A is given by 

A = H (5) 

and the corresponding scattering matrix is given hy S = K. 
Note that M is a Hermitian matrix. 

Definition 2: The linear quantum system (dJ is said to be lossless bounded real if the following 
conditions hold: 
i) F is a Hurwitz matrix; i.e., all of its eigenvalues have strictly negative real parts; 
ii) The transfer function matrix $(s) = H{sl — F)^^G + K satisfies ^(iuj^^iiuj) = I for all 

The following definition extends the standard linear systems notion of minimal realization to 
linear quantum systems of the form ([T]); see also lfT6ll . 

Definition 3: A linear quantum system of the form ^ is said to be minimal if the following 
conditions hold: 

i) Controllability. x^F = Xx^ for some A G C and x^G = implies x = 0; 



ii) Observability. Fx = \x for some A G C and Hx = implies a: = 0. 
The following theorem is an straightforward extension of Theorem 6.6 of [[T6l to allow for a 
general scattering matrix S. 

Theorem 2: A minimal linear quantum system of the form ([T]) is physically realizable if and 
only if the system is lossless bounded real. 

III. Singularly Perturbed Linear Quantum Systems 

A. General Singular Perturbations 

We now consider a class of quantum systems of the form ([T]) dependent on a parameter e > 
which are referred to as singularly perturbed quantum systems: 

dai{t) = Fnai{t)dt + Fi2a2{t)dt + Gidu{t)] 
da2(t) = -F2iai{t)dt -\ — F22a2(t)dt -\ — G2du(t); 
dy{t) = Hiai{t)dt + H2a2it)dt + Kdu{t). (6) 

This system can be re-written in the more standard singularly perturbed form (e.g., see ll26l ): 

dai{t) = Fiiai(t)dt + Fi2a2{t)dt + Gidu{t); 
eda2{t) = F2iai{t)dt + F22a2{t)dt + G2du{t); 
dy{t) = Hiai{t)dt + H2a2{t)dt + Kdu{t). (7) 

If the matrix F22 is non-singular, we can define the corresponding reduced dimension slow 
subsystem (e.g., see [|26ll ) by formally setting e = in (|7]) to obtain 

dai(t) = FQai(t)dt + Godu(t); 
dyit) = Hoai{t)dt + Kodu{t) (8) 

where 

Fq = Fii — F12F22 F21] 

Go = Gi — F12F22 G2', 

Hq = Hi — H2F22 F21] 

Ko = K-H2F22^G2. (9) 



This is the singular perturbation approximation to the system ©. We are interested in whether 
the reduced dimension quantum system ([8]), dll) is physically realizable if the singularly perturbed 
quantum system ([7]) is physically realizable for all e > 0. One approach to addressing this 
question might be to apply Theorem [21 and indeed, we can obtain the following theorem which 
is the main result of the paper: 

Theorem 3: If the singularly perturbed linear quantum system ([7]) is physically realizable for 
all e > and the matrix F22 is non-singular, then the corresponding reduced dimension quantum 
system ([8]), ^ is such that the matrix Fq has all of its eigenvalues in the closed left half of the 
complex plane and the transfer function matrix <l>o(s) = Ho{sI — Fo)~^Go + A'o satisfies 

%{iujy^o{iuj) = I (10) 

for all uj eR. 

The proof of this theorem is given in the appendix. 

Note that this result is not sufficient to prove the physical realizability of the reduced dimension 
quantum system ([8]), dH) since the application of Theorem [21 requires that the system realization be 
minimal and hence the conditions of Theorem [2l will only be satisfied if the matrix Fq is Hurwitz. 
However, the properties established in this theorem indicate that in all but pathological cases, 
the singular perturbation approximation will yield a physically realizable reduced dimension 
system. These pathological cases can be detected by testing the eigenvalues and minimality of 
the reduced order system. In addition, the following example shows one such pathological system 
in which the singular perturbation approximation is not strictly Hurwitz and not minimal. 
Example We consider a singularly perturbed quantum linear system of the form ([7]) where 
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For each e > 0, we calculate the characteristic polynomial of the matrix F^ 
to be 



-^11 -^12 

F2i/e F22/e 
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From this, it follows using the Routh-Hurwitz criterion that the matrix F^ is Hurwitz for all 

/ 



e > 0. Furthermore, it is straightforward to verify that the matrix 6^ 
the conditions 

F,e, + Q,F} + G,Gl = 0; 
G, + e,H^ = 



I/e 



> satisfies 



(11) 



where G^ 



G2/e 



and H 



Hi H2 



. Hence, it follows from Theorem \T\ that this 



singularly perturbed quantum system is physically realizable for all e > 0. Furthermore, it 

follows from ([TT]) that this system is in fact minimal for all e > 0. However, when we consider 

1 

which is not Hurwitz. 



-1 



the reduced order approximate system, we calculate Fq = 

Also, Go = 0, i^o = 0, Kq = —I and thus, the reduced order system is not minimal. 

This example shows that a stronger result than Theorem [3l which guarantees minimality 
and Hurwitzness of the approximate system, cannot be obtained in the general case. In the 
next subsection, we consider a special class of singular perturbations for which the physical 
realizability of the reduced dimension system can be guaranteed. 



B. A Special Class of Singular Perturbations 

We now consider a special class of singularly perturbed physically realizable quantum systems 
of the form ^ defined in terms of the matrices S, A and M in Definition [H Indeed, we consider 
the case in which = /, 



A 



Ai ^A2 



■ M 



j=Ml, IM22 



for all e > where S'^S = I, and Mu and M22 are Hermitian matrices. Then, substituting these 
values into Q, we obtain the following linear quantum system of the form ([T]): 

dai{t) = - ( -AJAi + iMn ) ai{t)dt 

-^ (^^1^2 + «Mi2 j a2{t)dt - A\Sdu{t); 

da2{t) = -^ (^^2Ai + iMl2J a^{t)dt 

— f-A^As + iM22) a2{t)dt - —A\Sdu{t)] 

dy{t) = kiai{t)dt+—=k2a2{t)dt + Sdu{t). (12) 

If we make the change of variables d2{t) = 4^02 (t), this leads to the following singularly 
perturbed quantum system of the form Q: 

dai{t) = - ( -AJAi + zA/ii I ai(t)dt 



-A|A2 + iMi2 ] d2{t)dt - A\Sdu{t); 
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- (-AlAi + iMl2]ai{t)dt 

- (2A2A2 + ^M22 j d2{t)dt - -k\Sdu{t)] 



da2(t) 



dy{t) = Aiai{t)dt + A2d2{t)dt + Sdu{t). (13) 

Note that even though we formally let e — )■ in the singular perturbation approximation, this 
state space transformation can be applied for each fixed e > 0. Then, for the singularly perturbed 
linear quantum system (fT3l) . we can obtain the corresponding reduced dimension approximate 
system according to equations ([8]), dH). 

The following result is obtained for singularly perturbed linear quantum systems of the form 
(fT3l) . This result can also be derived from the general nonlinear results presented in the papers 
[|23l . [|25l . However, this result for the linear case is included here for the sake of completeness. 

Theorem 4: Consider a singularly perturbed linear quantum system (fT3l) which is physically 
realizable for all e > and suppose that the matrix — ( IA2A2 + iM22 ) is nonsingular. Then the 
corresponding reduced dimension approximate system defined by equations ([8]), dH) is physically 
realizable. 



The proof of this theorem is given in the appendix. 



IV. Illustrative Example 

We consider an example from quantum optics involving the interconnection of two optical 
cavities as shown in Figure \T\ Each optical cavity consists of two partially reflective mirrors 
which are spaced at a specified distance to give a cavity resonant frequency which corresponds 
to the frequency of the driving laser; e.g., see ifTSl . EOll . In practice, optical isolators would 
also need to be included in the optical connections between the cavities to ensure that the light 
traveled only in one direction. 
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Fig. L A linear optical quantum system. 



Here Ki and K2 are the coupling parameters of the first cavity and 7 is the coupling parameter 
of the second cavity. These parameters are determined by the physical characteristics of each 
cavity including the mirror reflectivities. The QSDE of the form ([T]) describing this quantum 
system is as follows: 



d 



a 
a 



dyi 



-VK27 



a 
a 



dt 



/Ri + VK2 

/7 



du 



2, 



a 
d 



dt + duo. 



(14) 



We wish to consider the reduced dimension approximation to this system which is obtained by 
letting 7 — !■ 00. This corresponds to the case in which the mirrors in Cavity 2 are perfectly 



reflecting and so there is a direct optical feedback from the output ?/2 of Cavity 1 into the input 
ui of Cavity 1. If we let e = 4, it is straightforward to verify that this system is a system 
of the form ^ with Q = I, S = I, Ai = ^JK[ + v^, A2 = 1, Mn = 0, M22 = 0, 
Myi = f {VKI — y/K^). With the change of variables a = y/^a = -j^a, the system becomes 



a 
ea 






^ 



a 
a 



dt 



^ 



/^ 



du2; 



dyi 



^1 



^ 



a 
a 



dt + du2 



(15) 



which is a singularly perturbed quantum system of the form (|7]). Hence, the corresponding 
reduced dimension slow subsystem ([8]), dH) is given by 



da 
dyi 



fiyi+J^l ^ ^f^^\ ^^^ + (^^^ _ ^A^j 



rfWo 



A'2 ) adt — du2. 



(16) 



Since the system (fT5l) satisfies the conditions of Theorem IH it follows from this theorem that 
the system (fT6l) will be physically realizable. This can also be verified directly by noting that 
the system (fT6l) satisfies the conditions of Theorem [T] with = 1. 

Note that for this example, if A'l = K2, then the reduced dimension quantum system is 
uncontrollable, unobservable and has a pole at the origin. 



V. Conclusions 

In this paper, we have considered the physical realizability properties of the singular pertur- 
bation approximation to a class of singularly perturbed linear quantum systems. These results 
may be useful in the modeling of linear quantum systems such as gravity wave detectors where 
a simplified model is required without sacrificing physical realizability. 
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Appendix 

Proof of Theorem \3\ 

If the singularly perturbed quantum system ([7]) is physically realizable for all e > 0, then it 
follows from Theorem [U that for all e > 0, there exists a matrix > such that the matrices 



-^11 Fi2 

-F21 -F22 



G. 



',G2 



H 



Hi H2 



; K 



satisfy the conditions ([3]). Hence, it follows from the first of these equalities and Fact 12.21.3 of 
ll30l that matrix F^ has all of its eigenvalues in the closed left half of the complex plane for all 
e > 0. Then, using a standard result on singularly perturbed linear systems (e.g., see Theorem 
3.1 on page 57 of ll26l ') it follows that the matrix Fq has all of its eigenvalues in the closed left 
half of the complex plane. 

With the matrices F^, G^, H and K defined as above, it follows by a straightforward but 
tedious calculation that we can write the transfer function $e(s) = H(sl — F^)^^G^ + K in the 

form: 

-1 



^As) 



Hn + HoFm I 



'2-^22 



es 



'21 



sI-F. 



'12 



-^22 

es 



'21 



x\Go + F, 



12 



F22 
es 



G, 



^1 G2 

es 



+Ko + H2F,-,' il 

where the matrices Fq, Gq, Hq, Kq are defined as in Q. 

Now for small values of e > 0, we can approximate the term {l — ^) in the above 
expression as follows: 

I-^y' = -esF,~,' + Oie^). 

From this, it follows that we can write 

$,(s) = {Ho~esH2F22^F2i){sI-Fo) 
X (Go - esFi2F^,'G2) 
+K^-esH2F^^^G2 + 0{e^). 



-' [I + esFi2F2-2'F2i{sI - Fo)-'] ' 



From this and some further straightforward manipulations and simplifications, we can obtain 

$,(3) = <l>o{s)-es{Ho{sI-For'F^2 + H2F,f)F^^'{F2,{sI-Fo)-'Go + G2) 

+0{e'). (17) 

Now using the fact that the matrices F^, G^, H and K satisfy the conditions ([3]), we will show 
that transfer function matrix $e(s)is unitary at all frequencies. Indeed using ([3]), we have for all 

a; gM 

H {iul - F,y^ G,G\ {-tool - Fjy^ H^ 

= H {iooi - F,y^ [{tool -F,)e + e {-iui - fI)] {-iui - Fiy^ h^ 

= HQ {-luI - FlY^ H^ + H {iui - F,y^ QH^ 

= -KG\ {-iui - F^y^ H^ -H {iui - F,y^ G,KI 

Now using the third equation of ^ and the fact that K is square, we have for all a; G M 

= I-KK^ -KG\{-iuI-Fiy^ H^ -H{iuI-F,y^G,K^ 

-H {iui - F,y^ G,G\ {-lul - Fiy^ H^ 

= I - ^^{iu)^^{iuy. 

Therefore, since $e(iu;) is square we have 

^,{iuy^,{iu) = iyueR (18) 

for all e > 0. Hence, it follows from (fTTI) and the fact that (fTSi) holds for all e > that we must 
have 

^o{iuy^o{iu) = I 

for all a; G M. This completes the proof of the theorem. □ 

Proof of Theorem |?] 

For the singularly perturbed linear quantum system (fT3l) . it is straightforward but tedious to 
verify that the corresponding reduced dimension slow subsystem ([8]), ^ is given by 

dai{t) = - liM+-A^A]ai{t)dt-A^Sdu{t); 

dy{t) = Kai{t)dt + Sdu{t) (19) 



where 

-1 



A = Ai - A2 (^AJAs + ?M22 j (^A^Ai + tMlA ; 
S = 5 - A2 ('iAt2A2 + ^Mas") A^^; 
M = Mn + -AjAa (^^IM - ^Msa J M22 (^^l^^ + ^^22 ) A^Ai 

+M.(iA5A,-.M.)"M.(iAtA..«.)"M. 

-^Mi2 (^A^Aa - ^M22 j At2A2 (^A5A2 + 2M22 j A^Ai 

-^AlA2 (^A^2A2 - ^M22 J A5A2 ( ^A5A2 + ^M22 ] MiV (20) 

Furthermore, it is straightforward to verify that S'^S = I and the matrix M is Hermitian. Hence, 
it follows from Definition [T] that the system (fT9l) is physically realizable with the matrices A, S, 
M defined as above and with 6 = /. This completes the proof of the theorem. □ 
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